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Very Large Array Observations of Disturbed Ion Flow 
from the Plasmasphere to the Nighttime Ionosphere 

J. F. Helmboldt^ & H. T. Interna^ 

We present the results of a multi-scale analysis of TEC fluctuations using a roughly 
five-hour observation of the bright radio source Virgo A with the Very Large Array (VLA) 
at 74 MHz in its B configuration. Our analysis combines data sensitive to fine-scale 
structure (~10 km and < 0.001 TECU in amplitude) along the line of sight to Virgo A as 
well as larger structures (hundreds of km) observed using several (~ 30) moderately bright 
sources in the field of view. The observations span a time period from midnight to dawn 
local time during 1 March 2001. Several groups of magnetic eastward directed (MED), 
wavelike disturbances were identified and determined to be located within the 
plasmasphere (2.1 < L < 2.9). We have also detected evidence of non- wavelike structures 
associated with these disturbances which are propagating roughly toward magnetic north. 
These likely represent a non-uniform density flow from the plasmasphere toward the 
nighttime ionosphere. AE and Kp indices and GPS TEC data indicate that during the 
observations, there were low levels of geomagnetic activity accompanied by somewhat 
localized depletions in ionospheric density. Thus, the observed plasmaspheric disturbance 
may be part of a flow triggered by these ionospheric depletions, likely associated with 
forcing from the lower atmosphere which is typically more prominent during quiet 
geomagnc;tic conditions. In addition, we have also observed several roughly westward 
directed and southeast directed waves located within the ionosphere. They are coincident 
in time with the plasmaspheric disturbances and may be related to the deposition of 
material onto the nighttime ionosphere. 



1. Introduction 

Because of the negative impact on satellite-based 
communication, satellite operation, and GPS preci- 
sion, disturbances within both the ionosphere and 
plasmasphere as well as interactions between the 
two have generated much interest and investigations. 
These disturbances run the gamut from geomagnetic 
storms, substorms, streamers, ionospheric scintilla- 
tions, storm enhanced densities, and traveling iono- 
spheric disturbances (TIDs), to name a few. Con- 
sequently, detailed explorations of these phenomena 
have been conducted using a wide variety of instru- 
ments including satellites, magnetometers, ionoson- 
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des, radar arrays, and GPS and other satellite beacon 
receivers. 

Within this field of study, radio-frequency inter- 
ferometers are a relatively underused resource. De- 
signed to operate as interferometric telescopes, such 
arrays arc affected by the ionosphere and plasmas- 
phere in the same way as ground-base satellite re- 
ceivers. Consequently, much effort has been devoted 
to compensating for these effects, especially within 
the VHF regime where the impact is much greater. 
These corrections contain information about density 
fluctuations within the ionosphere and/or plasmas- 
phere, being sensitive to the gradient of the total 
electron content (TEC) rather than the TEC itself 
[see, e.g., Jacobson and Erickson, 1992a; Cohen and 
Rottgering, 2009; Helmboldt et al, 2012a]. When 
operating within the VHF regime and observing a 
bright source, interferometers can measure differ- 
ences in TEC, or STEC between pairs of antennas 
to a precision of nearly 10"^ TECU (1 TECU= 10^^ 
m~^) and TEC gradients to a precision of ~ 10~^ 
TECU km-i [see, e.g., Helmboldt et al, 2012a]. 
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Because of its relatively stable electronics, large 
collecting area, and available VHF system (bands 
at 74 and 330 MHz), the Very Large Array (VLA; 
located at 34° 04' 43.497" N, 107° 37' 05.819" W) 
has been used on a somewhat limited basis to study 
the ionosphere and plasmasphcrc. Seminal investiga- 
tions by Jacobson and Erickson [1992a] and Jacob- 
son and Erickson [1992b] used 330 MHz observations 
with the VLA to identify and explore several wave- 
like phenomena identified within the data. Follow-up 
work done by Hoogeveen and Jacobson [1997] showed 
that a new class of waves found by Jacobson and Er- 
ickson [1992a] were located within the plasmasphere. 
They were shown to essentially be density fluctu- 
ations formed perpendicular to magnetic field lines 
at Mcllwain L-parameters between about 1.8 and 4 
which co-rotated with the plasmasphere with a mod- 
erate amount of (mostly) westward convection. In 
addition, a statistical description of the typical, pre- 
dominantly turbulent behavior of the ionosphere over 
the VLA using a large 74 MHz survey [Cohen et ai, 
2007] was presented by Cohen and Rottgering [2009]. 
Recently, Helmboldt et al. [2012b] detailed new spec- 
tral analysis techniques that enabled the simultane- 
ous identification and characterization of wave-like 
disturbances within VLA VHF data from scales of a 
few kilometers up to more than 100 km. 

Here, we detail a new interferometer-based ap- 
proach that builds upon this previous work of Ja- 
cobson and Erickson [1992a], Hoogeveen and Jacob- 
son [1997], and Helmboldt et al. [2012b] by using a 
~ 5 hour 74 MHz observation of a very bright VHF 
source as well as several moderately bright sources 
spread over a 15° area on the sky (or, about 80 km 
at an altitude of 300 km). This new technique al- 
lows for a more thorough evaluation of the environ- 
ment of fluctuations on medium to large scales than 
can be obtained with observations of a single source. 
This multi-source method also has better precision 
than would be possible with a similar GPS-based 
approach. When combined with simultaneous data 
for a single bright source in the field of view, we 
demonstrate that a unique characterization of both 
plasmasphcric and ionospheric wave properties from 
scales of a few to hundreds of kilometers can be at- 
tained. In §2, we describe the observations and the 
processing/analysis of the data with more detailed 
descriptions included as Appendices. In §3, we detail 
the properties of plasmasphcric disturbances within 
the data. In §4, we similarly describe ionospheric 



disturbances contained within the data. We include 
additional contemporaneous AE index and GPS data 
while interpreting the results in §5 and summarize in 
§6. 

2. Data and Analysis 

2.1. Observations and Initial Calibration 

The data-set we have chosen to analyze (VLA pro- 
gram number AD441) is a roughly five hour long ob- 
servation of the bright source Virgo A, or "Vir A" 
(also commonly known as 3C274) with the VLA at 
74 MHz. Vir A is among the brightest cosmic VHF 
sources with a flux density at 74 MHz of about 1,700 
Jy (1 Jy = 10^2^ W m"2 Hz^^). The observations 
were taken in the post-midnight hours local time of 
1 March 2001 starting at 06:08 UT using the VLA 
"4-band" system [Kassim, et at, 2007]. While obser- 
vations with this system are generically referred to 
as 74 MHz observations, the central frequency was 
more precisely 73.8 MHz with a 1.5 MHz wide band- 
pass. In this case, dual polarization was also used. 
The Kp index over this time period was 1-2 and the 
solar activity was fairly average (^'10.7wl30 SFU). 
The VLA was in its B configuration with the anten- 
nas spanning an area about 11 km in diameter and 
laid out according to the diagram displayed in Fig. 
1. 

2.2. Spectral Analysis 
2.2.1. Self-calibration Data 

The calibration, processing, and spectral analysis 

of the observations arc described in detail in Appen- 
dices A-B. Briefly, we used techniques developed by 
Helmboldt et al. [2012a] to convert the observations 
of Vir A to measurements of the TEC gradient time 
series for each VLA antenna. This conversion relies 
on a technique known as "self-calibration" which is 
discussed in detail in Cornwell and Fomalont [1999]. 
Hereafter, use of the term self-calibration in refer- 
ence to any results presented implies that they were 
based on these measurements. Two Fourier-based 
approachs were used to convert the TEC gradient 
time series (temporal sampling of 6.59 seconds) into 
mean spatial frequency spectra to identify distur- 
bances that may be present within the data. The 
flrst technique is described in detail by Helmboldt et 
al. [2012b] which uses the phases of the temporal 
Fourier-transforms to determine a spatial frequency, 
^, and direction for each temporal mode. Because it 
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is based solely on phase information, we refer to this 
technique as "phase-based." Within the phase-based 
technique, the determined direction of each mode is 
that of the line along which the phases of the Fourier 
transforms of the two TEC gradient components (i.e., 
north-south and east-west) agree. This gives this ap- 
proach the ability to sense the direction of motion 
for non-wavelike disturbances (see Appendix B for a 
more detailed discussion). 

The second method uses a straightforward dis- 
crete Fourier transform (DFT) of the antenna-based 
temporal Fourier transforms to map the spectral 
power as a function of north-south and east-west 
spatial frequency. This "DFT-based" method also 
employs an algorithm designed to mitigate the ef- 
fects of the rather poorly behaved impulse response 
of the Y-shaped VLA within this procedure (again, 
see Appendix B for more discussion). Because of 
the low spectral resolution of the relatively com- 
pact B-configuration VLA (11-km across), this tech- 
nique still only allows one to reliably estimate a sin- 
gle spatial frequency and direction for each tempo- 
ral mode (see Appendix B). Since it is a more stan- 
dard spectral approach, the DFT-based technique is 
mostly sensitive to wavelike disturbances for which 
it measures wavefront orientations and not necessar- 
ily directions of motion. This makes the two tech- 
niques complementary when both wavelike and non- 
wavelike disturbances are present. 

We have applied both techniques to the Vir A self- 
calibration data to produce two-dimensional spec- 
tral maps averaged within one-half-hour windows as 
functions of spatial frequency, ^, and azimuth an- 
gle measured clockwise from north. To aid with the 
identification of small-scale, small amplitude distur- 
bances, we have normalized these spectra by dividing 
each ^ bin by the total power over all azimuth angles 
within that bin. For convenience, we will refer to 
these specially normalized spectra as "azimuth spec- 
tra." 

The azimuth spectra derived from the Vir A data 

are displayed in Fig. 2. The results for the phase- 
based technique are shown as greyscale images and 
the DFT-based results are plotted as color-coded 
contours. The DFT-based azimuth spectra show 
several prominent features directed roughly toward 
magnetic east (azimuth of 102° for the VLA) be- 
tween local times of 01:30 and 03:00. Their orienta- 
tion implies that these disturbances are aligned along 
magnetic field lines. They also have relatively large 



speeds, mostly spanning the range 200-400 m s . 
These properties imply that they are likely located 
within the plasmasphere and are the same magnetic 
eastward directed (MED) disturbances discovered by 
Jacobson and Erickson [1992b] with the VLA. These 
were described in more detail by Hoogeveen and Ja- 
cobson [1997] as relatively long plasmaspheric irregu- 
larities stretched along field lines. After 03:00, west- 
ward and southwest directed waves begin to domi- 
nate the DFT-based azimuth spectra. 

The phase-based azimuth spectra also show detec- 
tions of the same phenomena, but with notably dif- 
ferent characteristics. Specifically, the MED waves 
are broken into several groups and the estimated 
azimuths for each group seem to "drift" along lin- 
ear structures toward the north with decreasing spa- 
tial frequency. This may indicate that the MED 
waves are associated with a significant amount of 
non- wavelike structure which is moving in a different 
direction, possibly along the field lines toward mag- 
netic north. We will discuss these results and the 
properties of the MED waves in more detail in §3. 
In contrast, the phase-based and DFT-based spec- 
tra agree rather well for the westward and south- 
westward directed phenomena seen after 03:00 local 
time, indicating they are likely dominated by wave- 
like structures. 

2.2.2. Field-based Calibration Data 

In addition to Vir A, several other moderately 

bright sources were also visible with the 15°-widc 74 
MHz VLA field of view. In particular, we identified 
29 sources from the VLA Low Frequency Sky Survey 
[Cohen et at, 2007] which could be used to measure 
the TEC gradient as a function of time along the 
line of sight to each source, using a technique called 
"field-based calibration" [see Cotton et ai, 2004, and 
appendix A]. We will use the term field-based cali- 
bration to refer to any results that were based on 
this analysis. Since the sources used within this tech- 
nique were not nearly as bright as Vir A (between 2.5 
and 10 Jy versus 1,700 Jy for Vir A), each of them 
could not be used to perform self-calibration. This 
is because the relative contributions of other sources 
to the measured visibilities is significant for fainter 
sources. For a source as bright as Vir A, these con- 
tributions are negligible. However, for a moderately 
bright source, one can measure the effect of the iono- 
sphere along the line of sight by making images of 
the source on relatively short time scales and mea- 
suring how its observed position changes with time 
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[see Appendix A and Cotton et ai, 2004, for more 
details] . 

These position shifts yield measurements of the 
mean TEC gradient over the array toward a partic- 
ular source, implying that these gradients have been 
smoothed with a roughly 11-km wide circular ker- 
nel (i.e., the diameter of the VLA B configuration). 
In order to have a high enough signal-to-noise ratio 
within each image to detect the source, each image 
used one minute of data, implying that the data were 
effectively smoothed in time as well. This spatial and 
temporal smoothing reduced the sensitivity of these 
measurements to small-scale, small amplitude fluc- 
tuations. Specifically, the measured spectral power 
drops to one half the actual power at a spatial fre- 
quency of 0.04 km~^ (or, a wavelength of 25 km) 
and/or at a temporal frequency of 26.5 hr~^ (or, a 
period of 2.26 minutes). The spectral response drops 
to zero at a wavelength of 11 km and/or a period of 
1 minute. 

When projected onto ionospheric heights (here, we 
used 300 km; see Appendix A), the pierce-points for 
these sources spanned nearly 100 km. This is illus- 
trated in Fig. 3 where wc have plotted the projection 
of the VLA antennas onto a height of 300 km cen- 
tered on each of the 29 sources. With such a distribu- 
tion of pierce-points, one can perform a straightfor- 
ward spectral analysis by applying a Fourier trans- 
form in three-dimensions, one temporal and two spa- 
tial. While such power spectrum "cubes" will be 
less sensitive to small-scale structures, they can be 
used to separate different Fourier modes with better 
fidelity than either the phase-based or DFT-based 
techniques applies to the self-calibration data. This 
is because both of these methods rely on the conver- 
sion of temporal modes into spatial ones to enhance 
spectral resolution, and if there are several waves 
with the same temporal frequency but different spa- 
tial scales, only the strongest of the waves will be 
detected. Thus, these two methods are quite com- 
plementary. 

We have developed a procedure for using the 29 
moderately bright sources to produce power spec- 
trum cubes which is detailed in Appendix B. We 
used this method to yield spectral cubes within the 
same half-hour intervals as were used for the azimuth 
spectra displayed in Fig. 2 for the self- calibration 
data. To display the wealth of information contained 
within these cubes, we have binned them by tem- 
poral frequency into four "channels," with frequency 



ranges of 0-6, 6-13, 13-21, and 21-30 hr'^. We have 
plotted the mean power spectra within each half-hour 
interval for each channel in Fig. 4 7 as grcyscale con- 
tours. For each of these "channel maps," we have 
plotted the shape of the impulse response in the 
lower left corner as an ellipse that represents the full 
width at half maximum (FWHM) of the impulse re- 
sponse function (IRF). We have also displayed maps 
of spectral power derived from the DFT-based tech- 
nique applied to the self-calibration data as false- 
color images with the channel maps. Each channel 
map shows varying levels of significant (i.e., > the 
size of the IRF) structure that, in some cases, overlap 
with features seen within the self-calibration-based 
maps. These structures will be discussed further in 
§4. 

3. Plasmaspheric Disturbances 

As demonstrated in §2.2.2, our analysis of the Vir 

A-bascd self-calibration data has revealed a group of 
MED waves similar to those detected by Jacobson 
and Erickson [1992b] and established by Hoogeveen 
and Jacobson [1997] to be plasmaspheric irregular- 
ities. We show the properties of these MED waves 
in more detail in Fig. 8. For each of the three half- 
hour intervals between 01:30 and 03:00 local time, we 
have plotted in the left panels the spectral power ver- 
sus spatial frequency projected along magnetic east, 
Cm e , for temporal modes with azimuths between 90° 
and 180° and speeds < 2000 km hr"! (556 m s'^). 
The latter criterion was used to exclude modes which 
were essentially noise. Most of the MED waves have 
wavelengths between 25 and 100 km with a strong 
concentration near roughly 30 km. This is very con- 
sistent with the range in wavelengths measured for 
similar plasmaspheric structures by Jacobson et al. 
[1996]. For the first two half-hour intervals, there 
also appears to be a distinct group of relatively large 
waves (wavelength ~ 140 km) which have spectral 
power nearly an order of magnitude stronger than 
the bulk of the MED waves. 

For the same temporal modes, we have also plot- 
ted spectral power as a function of speed projected 
along magnetic east, vme^ in the middle panels of 
Fig. 8. As noted before, most of the speeds are rela- 
tively high, ranging from about 150-500 m s~^. This 
is entirely consistent with the trace speeds observed 
for MED waves detected previously with the VLA 
[Hoogeveen and Jacobson, 1997]. Assuming that the 
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vast majority of the measured speed of any tempo- 
ral mode is simply the projection of the co-rotation 
speed of the plasmasphere along the vector normal 
to the wavefronts (i.e., magnetic east), we can esti- 
mate the heights and L-shells of these detected MED 
waves. The heights and L-shells corresponding to ob- 
served speeds of 100, 300, and 500 m s~^ are labeled 
in the middle panels of Fig. 8. They show that most 
of the MED waves arc consistent with plasmasphcric 
heights between 2000 and 7000 km with a prominent 
peak at about 3300 km for the middle time inter- 
val. For the geomagnetic latitudes of the VLA line 
of sight at these heights and times, this peak height 
corresponds to L= 2.27. The median L-shell for the 
VLA-detected MED waves detailed in Hoogeveen and 
Jacobson [1997] was 2.2, quite consistent with this es- 
timate. However, we note that the separate group of 
large/strong waves seen within the first two time in- 
tervals have speeds <80 m s~^. These are therefore 
likely ionospheric disturbances whose exact heights 
are less certain due to the influence of other factors 
on their motion such as the wind within the F-region. 

As stated in §2.2.2, our data are entirely con- 
sistent with the model of Jacobson et al. [1996] 
and Hoogeveen and Jacobson [1997] for MED waves 
that describes them as plasmaspheric irregularities 
stretched out along field lines. To estimate the zonal 
angular sizes of these disturbances as they map to 
lower altitudes, we have used the estimated heights 
and spatial frequencies to compute angular spatial 
frequencies, ^ang; in the magnetic eastward direction. 
The power as a function of ^ang is plotted in the right 
panels of Fig. 8. From these plots, one can see that 
most of the structures occupy scales between about 
0.1° and 0.5° with noticeable peaks in the second 
and third half-hour intervals corresponding to scales 
of about 0.3° and 0.2°. 

As noted in §2.2.2, the results from the phase- 
based spectral analysis technique show evidence for 
non-wavelike structure associated with the MED 
waves. This is illustrated by the approximately lin- 
ear features seen near magnetic east in Fig. 2. We 
have postulated that this may indicate movement of 
non- wavelike structures northward of magnetic east. 
Since these structures are likely associated with the 
plasmaspheric MED waves, the most logical direc- 
tion of motion for such structures would be toward 
magnetic north along the field lines. 

To demonstrate that this provides a good descrip- 
tion of the phase-based results, we have constructed a 
"toy" model of these irregularities. This model con- 



sisted of 10^ Gaussian TEC fluctuations with ran- 
domly generated sizes with a uniform distribution 
ranging from 15 to 40 km (here, "size" is (T\/27r, 
where a is the standard deviation for the Gaussian). 
The amplitudes of the Gaussians were drawn ran- 
domly from a Raylcigh distribution with a mode of 
unity (arbitrary units). The amplitudes of these 
Gaussians were then modified by a cosine density 
perturbation with a wavelength of A = 30 km and 
an azimuth of 100° by multiplying by a factor of 
[1 -|- 0.5cos(27rrG ■ fyj/X)] where fo is the position 
vector of each Gaussian and is a unit vector per- 
pendicular to the cosine wavefronts. The initial field 
of fiuctuations is displayed in the left panel of Fig. 
9 with the positions of the VLA antennas plotted as 
points in the center for reference. 

Within the model, both the Gaussian fluctuations 
and the cosine plane wave were given co-rotation 
speeds of 100 m s~^ due east. The Gaussian fluc- 
tuations were also given an additional flow veloc- 
ity parallel to the cosine wavefronts of 200 m s"-'^. 
We then simulated VLA measurements by advanc- 
ing the field of fluctuations using these constant ve- 
locities with 10-second time steps over a one-hour 
time period. We ran the simulation twice, once with 
the flow velocity directed toward magnetic north, 
and once toward magnetic south. As with the ac- 
tual self-calibration data, we fit a second order two- 
dimensional polynomial to the simulated 5TEC mea- 
surements (see Appendix A) and applied both the 
phase-based and DFT-based spectral analysis meth- 
ods to the results. As with the real VLA data, we 
used a sliding half-hour wide boxcar, thus yielding 
dynamic spectra over the central half-hour of the 
simulated observations. 

The results for the phase-based and DFT-based 
analyses for both the northward and southward fiows 
are shown in the right panels of Fig. 9 as images of the 
spectral amplitude as a function of spatial frequency 
and azimuth. As one would expect, the DFT-based 
approach successfully recovers the properties of the 
cosine wave with little evidence of the Gaussian fluc- 
tuations, albeit with a somewhat overestimated size 
for the southward flow (see the bottom right panel 
of Fig. 9) . Conversely, the phase-based spectra show 
ample evidence of the influence of the Gaussian dis- 
turbances. For both the northward and southward 
flows, the measured azimuths skew toward the direc- 
tion of motion for the Gaussians toward lower fre- 
quencies where they tend to dominate the spectral 
power. They form roughly linear structures very sim- 
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ilar to those evident in Fig. 2, even reproducing the 
"overshoot" to more southeastward azimuths (i.e., 
> 100°) for higher spatial frequencies seen within 
some of the features in Fig. 2. 

The exact shapes of these spectral features de- 
pend on several factors such as the Gaussian fluc- 
tuations' distributions of sizes and amplitudes and 
their flow speed as well as the number, sizes, and 
strengths of the field- aligned cosine waves. It is un- 
clear to what degree the structures detected within 
our phase-based spectra constrain the properties of 
the flow(s) and associated MED waves. However, 
it is clear that our toy model has demonstrated 
that a substantial flow of non-uniform density mate- 
rial along field lines toward magnetic north provides 
a adequate explanation of the structures observed 
within the phase-based azimuth spectra. Given that 
the VLA lines of sight throughout the observations 
probed latitudes above the geomagnetic equator, this 
implies flow of material down flux tubes. This is 
consistent with the fact that the observations were 
conducted at night when some refilling of the iono- 
sphere from the plasmasphere is expected [see, e.g.. 
Park and Banks, 1974]. 

4. Ionospheric Waves 

In addition to the significant presence of plasma- 
pheric activity within the data, there is also evidence 
of ionospheric disturbances, chiefly among the field- 
based calibration data. Since the field-based data are 
based on observations of several sources with non- 
parallel lines of sight, it is much more effective at 
probing ionospheric rather than plasmaspheric struc- 
tures. For instance, at a height of 300 km, the mean 
separation among the lines of sight to the 29 sources 
used is about 12 km. This implies an approximate 
Nyquist sampling limit of a minimum of 24-km-sized 
structures. Indeed, the channel maps shown in Fig. 
4-7 show little if any structure beyond spatial fre- 
quencies of 1/12 = 0.042 km~^. If we increase the 
assumed height to 3000 km, the approximate height 
of the plasmaspheric MED waves described in the 
previous section, the minimum probable scale like- 
wise increases to 240 km. This is much too large to 
detect any of the plasmaspheric MED waves which 
have wavelengths < 100 km. 

The channel maps in Fig. 4-7 give a glimpse of the 
large amount of information about ionospheric dis- 
turbances contained within the field-based spectral 



cubes. At the beginning of the observing run, rel- 
atively low-frequency, northeastward-directed waves 
are prominent with weaker, higher-frequency waves 
directed toward the south/southeast. These all have 
wavelengths between 50 and 100 km. Starting at 
a local time of 01:47, the ionospheric MED waves 
described in §3 appear over a wide range of tempo- 
ral frequencies, as high as 20 hr~^. Unlike the self- 
calibration data, they are seen within the field-based 
data through the remainder of the observing run. 
Westward/southwestward directed waves are seen to 
overlap with those detected with the self-calibration 
data starting at a local time of about 03:30, most 
prominently for frequencies >6 hr~^. However, un- 
like the self-calibration data, roughly westward (some 
northwestward, some southwestward) waves are ap- 
parent at all times starting at the third time bin cen- 
tered at a local time of 01:47. They span nearly the 
entire range of temporal frequency and have wave- 
lengths of between 50 and 100 km, similar to the 
ionospheric MED waves. 

5. Interpretation 

The relatively low Kp indices for the time of our 
observations and the days preceding it (< 3) indicate 
that the plasmaspheric disturbances wc observed are 
not related in any way to geomagnetic storm activity. 
They instead appear to be the same plasmaspheric 
irregularities elongated along magnetic field lines dis- 
covered by Jacobson and Erickson [1992b] and de- 
scribed in more detail by Jacobson et al. [1996] and 
Hoogeveen and Jacobson [1997]. However, unlike the 
methods employed by these authors which rely on fit- 
ting a single plane wave to the VLA ^TEC measure- 
ments for a given time interval, our spectral analysis 
techniques can break the data up into many wave- 
like components. These techniques reveal that at a 
given time, there can be many plasmaspheric MED 
waves occupying a wide range of L-shells and sizes 
as well as ionospheric MED waves. In addition, our 
specialized phase-based spectral analysis method has 
revealed the existence of non-uniform density fiows 
of material down the flux tubes occupied by these 
irregularities toward the nighttime ionosphere. 

In order for such a flow to take place, there needs 
to be some depletion within the ionosphere. Of 
course, during the night, recombination naturally 
causes such a depletion. However, it is well estab- 
lished that during times of low geomagnetic activ- 
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ity, the density within the ionosphere can fluctuate 
substantially, between about 25% and 40% on scales 
of a few hours to several days [Forbes et at, 2000; 
Mikhailov et al., 2004, 2007]. It may be that during 
our VLA observations, larger than average nighttime 
depletions occurred as a result of quict-timc variabil- 
ity caused by forcing from the lower atmosphere via 
planetary waves, gravity waves, and/or tides [Chen, 
1992: Rishheth and Medillo, 2001; Liu et al, 2010]. 

In the upper panel of Fig. 10 we establish that in- 
deed, the level of geomagnetic activity was rather low 
during our observations. We have done this by plot- 
ting the AE index measure within 5-minute intervals 
on the date of the observation, 1 March 2001, and 
the mean within ztl5 days, both as functions of VLA 
local time. While evidence of substorm activity can 
be seen near sunset and midmorning, the AE index 
during the time of our observations was substantially 
lower that the typical value for that time of year (i.e., 
late winter/early spring 2001). 

To show that this low level of activity was accom- 
panied by ionospheric depletions, we obtained GPS- 
based vertical TEC maps for North America from 
the MIT Haystack Observatory Madrigal Database 
for 1 March 2001 and for 15 days before and after 
this date. We used these data to construct maps of 
vertical TEC at 5°x5° resolution for geomagnetic lon- 
gitudes between —125° and —75° and geomagnetic 
latitudes from 30° to 60° at 5-minute intervals. We 
then computed a 30-day mean map at each 5-minute 
interval and divided it into the corresponding map 
for 1 March 2001 to remove typical intraday vari- 
ations. We then computed the mean of the log of 
this ratio at each latitude and longitude during the 
time of our observations and have displayed the re- 
sult in the lower panel of Fig. 10. From this map, one 
can see that there were several persistent, localized 
depletions relative to what was typical for the late 
winter/early spring of 2001. There were also several 
locations of relatively enhanced TEC, but most of 
the area in the map shown seems to be below av- 
erage. This is qualitatively consistent with the idea 
that these variations are the result of forcing from 
the lower atmosphere which can have a non-uniform 
effect on the observed TEC [see, e.g., Liu et at, 2010]. 

For comparison, we have also plotted contours of 
spectral power as functions of estimated geomagnetic 
longitude and invariant latitude. A, for the plasmas- 
pheric MED waves described in §3 and Fig. 8. This 
shows that the L-shells occupied by these irregular- 



ities map well to relatively depleted regions of the 
ionosphere. It may then be that the irregularities 
themselves are tied to substructure within these de- 
pleted regions. Such structure may cause some flux 
tubes to be more depleted than others as they re- 
fill the ionosphere, causing wavclikc features as the 
ion density fluctuates as a function of geomagnetic 
longitude. If there are in fact significant structures 
on scales of a few tenths of a degree in geomagnetic 
longitude (see §3 and the right panels of Fig. 8), this 
could contribute to the formation of the observed 
plasmasphcric MED waves. The GPS-based data 
does not provide the level of spatial resolution re- 
quired to search for such structures. Our own VLA 
data has the ability to detect features on these scales, 
but if they exist, they are unfortunately overwhelmed 
by the plasmaspheric MED waves and ionospheric 
disturbances. Future studies based on satellite bea- 
con data may be able to better test this hypothesis. 

There are two classes of ionospheric disturbances 
coincident in time with the plasmaspheric distur- 
bances, field-aligned waves directed toward magnetic 
east and roughly westward directed waves. Both 
classes of waves have wavelengths between about 
50 and 100 km, and therefore can be categorized 
as medium scale traveling ionospheric disturbances 
(MSTIDs). However, they occupy a relatively large 
range in temporal frequencies, including higher fre- 
quencies than have been previously observed. With 
amplitudes ^0.01 TECU, they are also weaker than 
typical MSTIDs [see, e.g., Herndndez-Pajares et al., 
2006; Tsugawa et al., 2007). This is likely a reflection 
of the limitations in sensitivity and time resolution 
of the instruments used to characterize MSTIDs in 
past studies (typically GPS receivers). 

The fact that they occur at the same time as the 
observed plasmaspheric disturbances may indicate 
that the generation of these waves may be related 
to the deposition of material from the plasmasphcre 
onto the ionosphere. The MED, field- aligned waves 
may be generated in the topside ionosphere where 
other factors such as winds will have less of an ef- 
fect and the disturbances are more likely to remain 
aligned with the magnetic field. The westward di- 
rected waves may arise lower within the .F-region 
where the wind has more of an influence. According 
to the publicly available GPI data-driven runs of the 
TIEGCM code [Wang et al., 1999, currently avail- 
able at http://www.hao.ucar.edu/modeling/tgcm/], 
the zonal wind above about 350 km is generally west- 
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ward near the location of the VLA during the time 
of day and year of our observations. 

6. Conclusions 

Using the unique capabihtics of the VLA VHF 
system, we have developed a multi-scale approach 
to the spectral analysis of TEC gradients. Probing 
a wide range of spatial scales and amplitudes, we 
have identified both plasmaspheric and ionospheric 
waves present during a ~ 5 hour observation during 
the night of 1 March 2001. Contemporaneous data 
suggest that geomagnetic activity was relatively low 
during this time, implying that the plasmaspheric 
disturbances were not part of storm-related activity. 
Instead, combining the VLA data with contempora- 
neous AE index and GPS data, we have concluded 
that the field aligned waves in the plasmasphcre ap- 
pear to be associated with the flow of non-uniform 
density material along their wavefronts toward rel- 
atively depleted regions of the ionosphere. The low 
level of geomagnetic activity and the spatial distri- 
bution of ionosphere depletions hint that the deple- 
tions may be linked to density fluctuations caused 
by forcing from the lower atmosphere via planetary 
waves, gravity waves, and/or tides. We are currently 
conducting a more detailed investigation of this pos- 
sibility using seven years of GPS-based TEC maps, 
AE index data, and 74 MHz VLA data from the VLA 
Low- frequency Sky Survey [Cohen et ai, 2007]. 

We have also simultaneously detected relatively 
fast (roughly 150 350 m s~^) and weak (amplitudes 
<0.01 TECU) MSTIDs within the ionosphere. They 
were coincidence in time with the plasmaspheric dis- 
turbances and may have been generated by the de- 
position of material from the plasmasphcre onto the 
nighttime ionosphere. They were directed rough 
toward the west and toward magnetic east. This 
may be a reflection of the altitudes at which these 
waves were generated. The field-aligned waves may 
have been generated relatively high within the iono- 
sphere whereas the westward moving waves may have 
formed lower in the middle F-region where the wind 
has a larger effect. 

Appendix A: Data and Calibration 

Al. Initial Calibration 

The data were obtained from the National Ra- 
dio Astronomy Observatory (NRAO) data archive^ 
and initially processed within the Astronomical Im- 



age Processing System [AIPS; Bride and Greisen, 
1994, http://www.aips.nrao.edu]. Within AIPS, the 
complex bandpass response of each antenna was com- 
puted and the data were corrected for this. An initial 
calibration was performed using a model of the dis- 
tribution of the intensity of Vir A on the sky pre- 
sented by Kassim et al. [2007] (currently publicly 
available at http://lwa.nrl.navy.mil/tutorial/). This 
model was used to compute model values of the cor- 
relations of signals between pairs of antennas, or "vis- 
ibilities." This model visibility data-set was then di- 
vided into the observed visibilities and the antenna- 
based gains were solved for within each scan (i.e., a 
contiguous block of observing time) using standard 
AIPS tasks. This calibration technique is commonly 
referred to as "self-calibration"; its application and 
robustness are detailed in [Cornwell and Fomalont, 
1999]. In this instance, self-calibration was used over 
relatively long time intervals to correct for amplitude 
and phase variations in the responses of the individ- 
ual antennas which generally vary on relatively long 
time scales. This also has the effect of removing any 
slowly varying ionospheric contribution to the ob- 
served phases, hut this has little impact on our anal- 
ysis since we are focussing on variations on much 
shorter time scales (see below). 

A2. Self-calibration 

Following the initial calibration, a second round 

of self-calibration was performed, this time at much 
shorter time intervals of 6.59 seconds, the shortest 
possible for this data-set given the time-averaging 
used to measure the visibilities. The motivation for 
performing two separate rounds of self-calibration 
will become apparent in the following section where 
we will detail a separate calibration performed over 
the entire 74 MHz field of view. 

This second round of self-calibration was used to 
compute the difference in TEC between each an- 
tenna and the reference antenna [see Cornwell and 
Fomalont, 1999, and Fig. 1] according to Helmboldt 
et al. [2012a]. Unlike the data-set used by Helm- 
boldt et al. [2012a], we found no obvious jumps in 
the self-calibration determined antenna phases which 
would negatively impact the process of unwrapping 
the phases needed to compute the values of (5TEC as 
functions of time. Instead, we were able to simply un- 
wrap the phases as functions of time within each scan 
and polarization. The phases were then de-trended 
by subtracting smoothed versions of the phase time 
series using a 1/2 hour wide, sliding boxcar filter. 
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Near the edges of each scan, a linear fit to the phases 
as functions of time were also subtracted to miti- 
gate edge effects within the smoothing process. The 
boxcar was chosen to have the narrowest width that 
would not substantially dampen the longest period 
waves expected within the spectral analysis to be 
perform (see §4). We expect to see waves with peri- 
ods up to about one hour [e.g., Hemdndez-Pajares et 
al, 2006; Helmholdt et al, 2012b]. The subtraction 
of a version of the data smoothed with a 1/2 hour 
wide boxcar will reduce the amplitudes of such waves 
by a factor of about 2.75 which is acceptable given 
that these waves also tend to have relatively large 
amplitudes. 

After this, the two polarizations were averaged and 
divided by a factor of 114 [see, e.g., equation (3) 
of Helmboldt et al., 2012a] to put them in units of 
TECU. The rms deviations between the two polar- 
izations within bins of 10 time steps were used to es- 
timate the la uncertainties in the JTEC values. Fig. 
11 shows examples of (5TEC time series for antennas 
in the northern arm of the VLA (see Fig. 1). The 
typical error in the JTEC values is about 3 x 10~^ 
TECU. This is very similar to what was achieved by 
Helmboldt et al. [2012a] using a source that is about 
ten times brighter than Vir A with similar time aver- 
aging (6.67 seconds). This implies that the accuracy 
attained is largely influenced by the time averaging 
used in the observations and that there are real iono- 
spheric fluctuations on the order of 10~^ TECU on 
time scales of about 6.5 seconds. 

A3. Field-based Calibration 

To obtain some measure of ionospheric fluctua- 
tions over the entire 15° field of view of the 74 MHz 
VLA, we have also performed a procedure referred 
to as "field-based" calibration. The principles behind 
this method and its common application are detailed 
by Cotton et al. [2004] and Cohen et al. [2007]. In 
short, the first order effect of the ionosphere on an 
image of a cosmic source made from interferometric 
data is to shift its position. The observed position 
shifts of a particular source are proportional to the 
components of the TEC gradient along the line of 
sight to the source averaged over the span of the ar- 
ray and the time interval used to make the image. 
At 74 MHz, images of sources with flux densities 
of more than about 2.5 Jy can typically be made 
with 1-2 minutes of VLA data to the quality needed 
to accurately measure the positions of the sources. 



Field-based calibration uses a predetermined catalog 
of sources to choose a set of bright "calibrators" and 
makes images of each calibrator within 1 2 minute 
intervals, commonly referred to as "snapshots". The 
position offset from its expected location is then mea- 
sured for each source and a Zernike polynomial is fit 
to all of the position offset data for a particular time 
interval. The Zernike polynomials are then used to 
correct the entire field of view for ionospheric effects. 

For the observation of the field surrounding Vir 
A, an additional step needed to be performed be- 
fore running field-based calibration. Since Vir A is 
such a bright source, the dominant source of noise 
within an image made of any source in its vicinity 
is secondary peaks in the impulse response function 
(IRF), or "sidelobes" from Vir A. It was therefore 
necessary to subtract the contribution of Vir A from 
the observed visibilities before field-based calibration 
was implemented. This was done by using the second 
round of self-calibration to add the effects of the iono- 
sphere to the model visibility data-set computed with 
the model of the sky brightness distribution of Vir 
A. This altered model data-set was then subtracted 
from the visibility data. Since the self-calibration 
process was done in two steps, this subtracted data 
was still corrected for instrumental effects while the 
short time-scale ionospheric variations remained. 

Following the subtraction of Vir A, field-based cal- 
ibration was performed using the software package 
Obit (http:/ /www. cv. nrao.edu/~bcotton/Obit. html) 
using one minute snapshots of all calibrators within 
the field of view brighter than 2.5 Jy at 74 MHz, 
48 sources in all (excluding Vir A). The input 
source catalog was that produced by Cohen et al. 
[2007] from the VLA Low Frequency Sky Survey 
(VLSS) data. This run of field-based calibration 
was somewhat compromised by the presence of radio 
frequency interference (RFI) which is usually seen 
within 74 MHz VLA data, especially on shorter base- 
lines. However, a relatively new technique to sub- 
tract the effects of stationary sources of RFI from 
interferometric data has been developed by Athreya 
[2009] and implemented within Obit. This technique 
requires that relatively bright sources be subtracted 
from the data before it is used, otherwise it may 
subtract potions of the intensity of such sources at 
times when their visibility fringe-rates are similar 
to those of the RFI. Thus, we used the first run 
of field-based calibration to subtract out all sources 
detected within the field of view (this is automati- 
cally done within Obit) and ran the RFI subtraction 
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algorithm on the resulting residual data. We then 
ran an additional RFI flagging routine within Obit 
which flags data with unusually large amplitudes, es- 
pecially those with large Stokes V signals as few if 
any cosmic VHF sources have non-zero Stokes V in- 
tensities. Following this, the previously subtracted 
sources were added back to the data, and field-based 
calibration was performed a second time. On aver- 
age, the RFI subtraction and flagging reduced the 
noise in the calibrator snapshot images by 20%. 

Following the completion of field-based calibra- 
tion, we used the measured calibrator position shifts 
to compute the TEC gradient as a function of local 
time for each source. We limited ourselves to only 
those sources that were detected in their snapshot 
images more than half of the time so that the the 
resulting time series would be reasonably complete 
which will be important for subsequent spectral anal- 
ysis (see Appendix B) . There were 29 calibrators that 
met this criterion. The relative positions of their as- 
sociated ionospheric pierce-points are plotted in Fig. 
3 for a height of 300 km. Since we do not know a pri- 
ori what height most of the observed fiuctuations will 
occur at and the disturbances will likely span a range 
of heights, the choice of altitude is somewhat arbi- 
trary. We have chosen 300 km as an approximate 
location within the lower F region where medium 
scale disturbances are likely to be found. In Fig. 
3, the positions of the VLA antennas are also plot- 
ted around each calibrator position at a height of 300 
km to represent the amount of smoothing that has 
been eflFectively applied by using the position offsets 
to measure the local TEC gradients. 

After the position shifts were converted TEC gra- 
dients assuming a height of 300 km, they were de- 
trended in a similar manner to the self-calibration- 
derived (5TEC values for consistency. In particular, 
for each component of the gradient and at each time 
interval, a line was fit to the data for all time inter- 
vals within ±1/4 hour. Any data with absolute devi- 
ations from the fitted line more than three times the 
median absolute deviation (MAD) were then iden- 
tified and the fit was redone without them. If the 
absolute difference between the value at the current 
time interval and the final linear fit was <5xMAD, 
the fitted value was subtracted from it. Otherwise, 
this time interval was flagged as being possibly spu- 
rious. 

As an example of the results of this processing, 
we have plotted the de-trended north-south compo- 



nents of the TEC gradients as functions of time for 
all 29 calibrators in Fig. 12 with the source peak in- 
tensities printed in each panel. Note that these are 
the observed peak intensities which are lower than 
the sources' cataloged intensities because of the re- 
sponse of the individual antennas which decreases 
significantly with angular separation from the center 
of the field of view. Flagged time intervals are high- 
lighted in red. These data include error-bars com- 
puted using the errors in the position offsets. These 
were computed at each time interval according to 
Condon [1997], where the error in the source posi- 
tion is given by 

CTp « ^ae/Sp (Al) 

where a is the noise in the image. 6 is the angular di- 
ameter of the source, and Sp is the peak intensity on 
the image. Here, we have assumed that the sources 
arc unresolved and consequently have the same angu- 
lar diameter as the IRE, in this case, 76 arcseconds. 
The dominant source of noise within each snapshot is 
sidclobes from other sources within the filed of view. 
To estimate the typical value of a for the snapshots, 
we made an image of the full filed of view using all 
270 minutes of data with no deconvolution applied 
(i.e., without removing source sidelobes). The rms 
noise for this image was 0.042 Jy beam~^, implying 
a « 0.7 Jy beam~^ for each one-minute snapshot. 
Given the intensities of the calibrators used (see Fig. 
12), this further implies a typical TEC gradient un- 
certainty of about 10"^ TECU km'^ 

Appendix B: Analysis Methods 

We have performed separate spectral analyses on 
both the self-calibration and field-based calibration 
data. As we will detail below, the methods used are 
quite complementary and we include some compari- 
son of the two in the main body of this article. 

Bl. Self-calibration Data 

The spectral analysis of the de-trended 6TEC data 
derived from self-calibration using Vir A (see Fig. 
11 for examples) was performed using two comple- 
mentary techniques. The first follows that devel- 
oped by Helmholdt et al. [2012b] and the second uses 
a straightforward Fourier-style analysis. Both tech- 
niques start with a second order, two-dimensional 
polynomial fit over the array to obtain approximate 
values for the TEC gradient at each antenna's as- 
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sociated pierce-point. The locations of these pierce- 
points were computed assuming a nominal height of 
300 km. However, since the lines of sight from the 
antennas to the sources are essentially parallel, the 
effect of the assumed height on these computed po- 
sitions is marginal unless the observed source is at 
a relatively low elevation. During the VLA observa- 
tions, Vir A was no fewer than 45° above the hori- 
zon, indicating that the projected antenna separa- 
tions changed very little throughout the observations 
Helmboldt et al. [see 2012a] and had only a weak de- 
pendence on assumed height. For instance, the rela- 
tive pierce-point positions among the antennas for a 
height of 300 km differ by at most a few percent from 
those computed for a height of 3000 km. The time 
series for the coefficients of the polynomials fit using 
these pierce-point positions were then Fourier trans- 
formed (here, within a 1/2 hour sliding window) to 
obtain a complex spectrum of the gradient for each 
antenna. 

Bl.l. Phase-based Technique 

The two spectral analysis techniques both use the 
temporal Fourier transforms of the antenna-based 
TEC gradient time scries. The first technique from 
Helmboldt et al. [2012b] uses the phase information 
from the Fourier transform at each time step to esti- 
mate the weighted mean direction and magnitude of 
the wavenumber vector, k, for each temporal mode. 
These results were combined with the temporal fre- 
quencies to also estimate pattern speeds. This tech- 
nique uses both components (i.e., north-south and 
east- west) components of the TEC gradient. The 
technique assumes that the line along which the 
phases of the Fourier transforms of the two compo- 
nents agrees is the direction for that temporal mode. 
This is done by first fitting a plane to the Fourier 
transform phase of each TEC gradient component 
after unwrapping the phases along each VLA arm. 
Then, the parameters of the two planar fits are used 
to determine the line defined by the intersection of 
the two planes. 

It was shown by Helmboldt et al. [2012b] that this 
phase-based approach can mitigate the influence of 
wavefront distortions on the analysis of relatively 
large waves (i.e., those substantially larger than the 
array). However, it also has the potential to sense 
the direction of motion for non- wavelike patterns. 
Imagine, for instance, a distortion of arbitrary shape 
passing over the array. If one chooses a series of 
points along its direction of motion, the time series 



will be identical for each point except that they will 
be out of phase with one another. If the velocity is 
constant, then the phase difference between two of 
these points will just be proportional to the sepa- 
ration between them. This holds true for the TEC 
times scries as well as for the time scries of the two 
components of the gradient. Therefore, the "phase- 
based" technique of Helmboldt et al. [2012b] is capa- 
ble of estimating the direction of motion, rather than 
orientation, when non-wavelike disturbances are en- 
countered. 

B1.2. DFT-based Technique 

The second technique uses a two-dimensional, spa- 
tial discrete Fourier transform (DFT) of each of 
the TEC gradient components measured at each 
antenna^. Since the array is "Y" -shaped, one can 
imagine that the IRF for this DFT is rather poorly 
behaved. We show in Fig. 13 that this is in fact the 
case. The left panel of Fig. 13 displays the IRF for 
the first half-hour of VLA data. As one might ex- 
pect, there are substantial sidelobes in the regions 
between VLA arms extending to relatively large spa- 
tial frequencies. To mitigate the effect of these side- 
lobes, we have implemented a version of the clean 
algorithm [Cornwell et al, 1999]. This algorithm is 
often used with images produced using radio interfer- 
ometers to remove the effects of sidelobes associated 
with bright sources. Briefly, it is an iterative algo- 
rithm that identifies a peak in the image plane at 
each iteration and subtracts a scaled version of the 
IRF centered on that peak where the scale factor is 
usually the intensity multiplied by a predetermined 
gain value. The iterations are continued \int\\ a max- 
imum iteration limit is reached or some convergence 
criteria are met. Often, this is done interactively us- 
ing boxes drawn by eye to limit the locations of the 
model delta functions, or "clean components" to spe- 
cific regions. After the cleaning is finished, the clean 
components arc convolved with a more well-behaved 
IRF, usually a Gaussian whose parameters are de- 
termined from the main peak of the IRF. Following 
this, the re-convolved components are added back to 
the image. 

Our implementation is only slightly different since 

the spatial spectra we arc "cleaning" arc complex and 
not real valued functions like the images made with 
interferometers. Consequently, at each iteration, we 
used the spectral power to determine the location of 
the next clean component to be subtracted, and the 
clean components all had complex amplitudes. For 
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each clean component, a scaled version of the IRF 
was used that was computed using the same grid 
spacing as the spectra but with 2N + 1 pixels where 
A'' is the number used for the actual spectra (here, 
N = 256; see Fig. 13). This was done so that the 
IRF could be easily shifted to be centered directly 
on a particular spectral pixel. After trying several 
convergence criteria, maximum iteration limits, and 
gain values, we found that satisfactory results were 
obtained by simply running 20 clean iterations with 
a gain of 0.5 and no limits on the locations of the 
clean components (i.e., no "boxing"). This clean- 
ing was run separately on each temporal frequency 
and TEC gradient component within each time bin. 
After restoring the rc-convolvcd clean components to 
their respective spectra, the power spectra from both 
components of the TEC gradient were combined into 
a single gradient power spectrum for each temporal 
frequency and time bin. 

Fig. 13 shows examples of the results of the 
clean algorithm applied to the first half-hour of self- 
calibration data. In the right panels, we show the 
spectral power for the strongest temporal frequency, 
6.375 hr""'^, and for a randomly chosen frequency, 
75.374 hr~^, both before and after cleaning. From 
these, on can see that nearly all the structure seen 
within the raw spectral maps was from the IRF. It 
also shows that the overarching assumption of Helm- 
boldt et al. [2012b] that each temporal mode is well 
approximated by a single wavenumber vector is es- 
sentially correct. Because of this, and so we can 
more directly compare the two methods, we used the 
cleaned spectral power to computed the peak power 
and a single weighted (with the spectral power) mean 
value of the wavenumber vector for each temporal 
mode for each time step. Because it is a more con- 
ventional Fourier-based analysis, this "DFT-based" 
method is most sensitive to wavelike structures and 
can only measure the wavefront orientation, which 
dow not necessarily give the direction of motion. The 
importance of this distinction is made apparent in 
the discussion in §3. 

Because both techniques estimate a single 
wavenumber vector for each temporal mode, they es- 
sentially use both the apparent motion of Vir A and 
the motion of the observed disturbances to convert 
the properties of temporal oscillations observed with 
a compact array into spatial information. This al- 
lows a much higher spatial frequency resolution to be 
achieved than could normally be obtained with such 



a small array. For instance, the temporal spectral 
resolution, Av, is related to its spatial counterpart, 
A^, simply via the pattern speed, i.e., — Av/v. 
For the 1/2 hour wide window we have used here, 
Av « 2 hr~^. Therefore, for a disturbance traveling 
100 m s~^, A^ w 0.0055 km"'^ which is the equiva- 
lent of what can be obtained with an array with a 
diameter of nearly 200 km. The results from both 
techniques are discussed in more detail in §2. 

B2. Field-based Calibration Data 

While the self-calibration-bascd spectral analysis 
supplies results with remarkable sensitivity and spa- 
tial frequency resolution, it suffers from two signif- 
icant drawbacks. First, even though it uses a rela- 
tively long temporal baseline to achieve the equiva- 
lent of a long spatial baseline, the measurements are 
still confined to a narrow strip, in this case, about 
11 km wide. This makes it difficult to accurately 
estimate the direction and speed of relatively large 
disturbances in the presence of significant wavefront 
distortions. Second, wave properties are derived from 
the data for each temporal frequency. This implies 
that if there are several wave patterns with differ- 
ent sizes but similar periods, there is no way within 
this method to separate them from one another and 
all that can be estimated is a combined speed and 
direction for all of them. For isotropic phenomena 
such as turbulence, this may be acceptable, but for 
characterizing wave properties, it can be problem- 
atic, especially at temporal frequencies where there 
is no clearly dominant wave pattern. 

What is needed to supplement this analysis is an- 
other array that provides TEC gradient measure- 
ments over similar time intervals but spread out over 
a much larger area. The position offset data used in 
the field-based calibration method detailed in Ap- 
pendix A provides just such an array. As Fig. 3 
demonstrates, the distribution of calibrator sources 
on the sky provides and array of TEC gradient mea- 
surements that spans an area with a diameter of 
about 100 km or more (at a height of 300 km) with 
a quasi-random distribution of ionospheric pierce- 
points. With this data, it is therefore possible to 
perform a relatively straightforward Fourier analysis 
on the data without utilizing any "tricks" to retrieve 
spatial information about the observed TEC fluctu- 
ations. One may simply perform a Fourier transform 
in three dimensions, one temporal and two spatial, 
over a given time interval. This will yield spectral 
"cubes" so that different spatial modes with simi- 
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lar temporal frequencies may be separated from one 
another. We note that the field-based data has ef- 
fectively been smoothed in time (with a one minute 
window) and space (with an 11 km diameter circular 
window) so that higher frequency modes are damped 
out. However, the self-calibration data provide am- 
ple information about such modes, making the two 
analyses quite complementary. 
B2.1. Method 

To begin the analysis of the field-based data, we 
first put the data into 8 temporal bins. To minimize 
sidelobes in the temporal components of the Fourier 
transforms, we have weighted the TEC gradient time 
series within each bin by a unit amplitude Gaussian 
with a full width at half maximum (FWHM) of 1/2 
hour. We therefore constructed the temporal bins to 
be one hour wide, but separated by 1/2 hour so that 
all time intervals would have Gaussian weights > 0.5 
in at least one time bin. For each source, there were 
also time intervals where either the source's position 
could not be measured, or the data was otherwise 
flagged as spurious (see Appendix A). Such gaps in 
temporal coverage can also lead to increased side- 
lobes in frequency space. Wc therefore reduced the 
Gaussian weight of any time interval adjacent to such 
a gap by a factor of two to introduce some amount of 
tapering near these gaps to mitigate their effects on 
the impulse response. Using this weighting scheme, 
we performed a discrete Fourier transform (DFT) of 
each component of the TEC gradient time scries for 
each source within each of the 8 time bins at 60 fre- 
quencies up to the Nyquist limit of 30 hr~^. In each 
instance, wc also computed the amplitude of the im- 
pulse response at u = and normalized each complex 
temporal spectrum by this value. 

Following the generation of temporal spectra for 
the calibrators, the mean north-south and east-west 
position of each calibrator's pierce-point relative to 
that of the center of the filed of view at a height of 
300 km was computed for each time bin. These were 
then used with a DFT to compute a complex spec- 
trum on a grid of north-south and cast-west spatial 
frequencies for each time bin and temporal frequency 
and for each component of the TEC gradient. Again, 
these spectra were normalized by the amplitude of 
the IRF at a frequency of zero. Given the relatively 
sparse sampling provided by the distribution of cal- 
ibrators, each resulting spatial spectrum was signifi- 
cantly affected by sidelobes. We sought to mitigate 
these effects using an implementation of the clean 



algorithm as was done with the DFT-based method 
that was applied to the self-calibration data (see §B1 
and Fig. 13). As before, 20 iterations were used per 
temporal frequency and gradient component with a 
gain of 0.5. In Fig. 14, we show the mean gradi- 
ent power over all temporal frequencies at each time 
bin both before and after cleaning. For reference, we 
have also included images of the IRF power spectra. 
One can see that the cleaning process has substan- 
tially improved the fidelity of the spectral images. 

Unlike the self-calibration data, the multiple lines 
of sight used in this analysis are far from parallel. 
This means that each of them "views" any distur- 
bance from a different angle causing the observed 
amplitude of such a disturbance to vary across the 
VLA field of view. For instance, if there is an op- 
timum viewing angle for a particular disturbance, 
the detected amplitude for a given line of sight is 
(X cos(9ios), where 9ios is the angle between the ac- 
tual line of sight and the optimum one (i.e., the dot 
product between the unit vectors for the actual and 
optimum lines of sight). The effect of this on our 
analysis is that this position dependent weighting of 
the measured gradients alters the IRF from its ex- 
pected shape. To illustrate the magnitude of this 
effect, we show in the panels of Fig. 15 the observed 
IRF for values of dios ranging from 10° 80°. For this 
demonstration, we assumed the optimum view an- 
gle was toward the north (i.e., its unit vector had 
no east-west component). As expected, the most no- 
table change is for the largest 6io8 = 80°. The lower 
panels of Fig. 15 show this better where they dis- 
play the difference between the observed IRF (scaled 
to a maximum of unity) and the expected IRF. The 
maximum difference for Oiog = 80° is nearly 10%, 
but ranges from <1% to a few percent for the other 
values of 0ios- 

Therefore, when we clean the field-base spectral 
cubes, we are basically using the wrong IRF, leading 
to sidelobes that are not properly subtracted and 
even adding additional but weak sidelobes. Since 
there are typically several detected features within 
each field-based spectral cube which likely have a dis- 
tribution of optimum viewing angles, the net effect 
of this on our analysis is to increase the noise of the 
spectral maps. However, as Fig. 15 indicates, this is 
likely a relatively small effect which is typically on 
the order of 1% and only as high as 10% in extreme 
cases. 
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Notes 

1. https://axchive.nrao.edu 

2. Formally, an inverse DFT was used since a particular oscil- 
lating mode ~exp[27ri(i^t — x^x — 
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Figure 1. The layout of the VLA antennas during the 
observations. The reference antenna (see Appendix A) is 
highlighted in white. 
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Figure 2. For the Vir A self-calibration data, the mean 
power within bins of azimuth angle and spatial frequency 
with 1/2 hour bins for the phase-based (image) and DFT- 
based (contours) techniques (see §2 and Appendix B). An 
azimuth of 0° corresponds to north, 90° to east, —90° to 
west, and ±180° to south. Within each spatial frequency 
bin, the values have been normalized by the total over all 
azimuth angle bins to make high-frequency features more 
apparent. The local time at the center of each time bin 
is given in each panel. 
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Figure 3. The distribution of ionospheric pierce-points 
at an altitude of 300 km for all 29 calibrator sources used 
for field-based calibration (see §2 and Appendix A) rel- 
ative to that of Vir A (in red) when it was at its peak 
elevation in the sky. The individual VLA antenna po- 
sitions are plotted relative to each source's pierce-point 
with a circle drawn to indicate the area covered by the 
VLA. 
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Figure 4. Within 1/2 iiour windows, tiie mean spec- 
tral power as a function of spatial frequency for the 
field-based calibration data (see §2 and Appendix B) for 
temporal frequencies between and 6 hr~^ plotted as 
greyscale contours. The shape of the impulse response 
function (IRE) for each map is plotted in the lower left 
of each panel as an ellipse representing the full width 
at half maximum (FHWM) of the IRF. The underlying 
image within each panel is a map of the log of the total 
power from the self-calibration data using the DFT-based 
spectral analysis technique (see §2, Appendix B, and Fig. 
2). 
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Figure 5. The same as Fig. 4, but for temporal frequencies between 6 and 13 iir ^. 
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Figure 6. The same as Fig. 4, but for temporal frequencies between 13 and 21 iir ^. 
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Figure 7. The same as Fig. 4, but for temporal frequencies between 21 and 30 hr ^. 
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Figure 8. For the three half-our time bins where they 
are most prominent (see Fig. 2), the properties of the 
field-aligned plasmaspheric disturbances. The properties 
of each temporal mode from the DFT-based analysis with 
azimuths between 90° and 180° and speeds <2000 km 
hr~^ are plotted as separate points with magenta con- 
tours indicating the density of points within each panel. 
The left panel shows the spectral power as a function 
of the spatial frequency vector projected along magnetic 
east, ^ME- The middle panel plots the spectral power 
versus the velocity projected along magnetic east, vme, 
with corresponding estimates of the height, z, and L-shell 
computed assuming all of the motion is from co-rotation. 
The right panels combine the results from the other two 
columns of plots to show the spectral power as a function 
of angular spectral frequency, ^ong (i-e., the disturbance 
scale in geomagnetic longitude). 
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Figure 9. Results from the "toy" model discussed in 
§3. The left panel shows the initial distribution of 10* 
Gaussian fluctuations modified by a 30-km-wavelength 
plane wave aligned along the magnetic field. The posi- 
tions of the VLA antennas are plotted as magenta points. 
Results for both the phase-based and DFT-based anal- 
ysis techniques are shown in the right panels for fiows 
directed toward magnetic north (left panels) and toward 
magnetic south (right panels). In all four panels, the 
spectral amplitude is displayed as a function of spatial 
frequency (abscissa) and azimuth (ordinate). White dot- 
ted lines in all four panels indicate the actual azimuth 
and wavelength of the field-aligned plane wave. 
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Figure 10. Upper: The AE index averaged within 
5-minute bins for the date of the VLA observations, 1 
March 2001 (black curve), and averaged with ±15 days 
of the observations (grey curve) as functions of VLA local 
time. Lower: For the time range during the VLA obser- 
vations (06<UTcll), the mean of the log of the ratio of 
the vertical TEC on 1 March 2001 to the 30-day average 
(centered on 1 March 2001) within 5° x 5° bins in geo- 
magnetic latitude and longitude from the MIT Madrigal 
database. White contours represent the spectral power of 
field-aligned plasmaspheric disturbances as a function of 
geomagnetic longitude and estimated invariant latitude, 
A. 
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Figure 11. For antennas in the northern arm of the 
VLA, the de-trended values for 5TEC relative to the ref- 
erence antenna (see Fig. 1 and Appendix A) as functions 
of local time. The baseline lengths and median la un- 
certainties estimated using both polarizations are printed 
for reference within the panels. 
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Figure 12. For the 29 calibration sources, the de- 
trended and flagged (red points) north-south components 
of the TEC gradient (north is positive) as functions of 
local time. The peak intensity of each source is given 
in each panel in units of Jansky (1 Jy = 10~^^ W 
Hz-'). 
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Figure 13. An example of the DFT-based spectral anal- 
ysis on the first half-hour of self-calibration data. The 
left panel shows the impulse response from the DFT of 
the antenna positions. The panels to the right show the 
spectral power as functions of north-south and east-west 
spatial frequency for the brightest temporal frequency, 
6.375 hr-\ and one chosen at random, 75.375 hr ^ , both 
before and after the clean algorithm was applied. 
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Figure 14. Within 1/2 liour windows, tiie mean (over 
all temporal frequencies) TEC gradient power spectra de- 
termined from the field-based calibration data (see Ap- 
pendix A and Fig. 12 for examples) before and after the 
application of the clean de-convolution algorithm (see 
Appendix B). The IRF for each spectrum is also dis- 
played for reference. 
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Figure 15. Examples of the impact of different viewing 
angles of plane wave-like disturbances on the measured 
impulse response for the spectral analysis of field-based 
data. The upper panels show the measured IRF for dif- 
ferent values of the angle between the line of sight to 
the center of the field of view and the optimum viewing 
angle for a given disturbance, Oioa- The bottom panels 
show the difference between this IRF and the unaltered 
version (see §B2.1 for more details). 



